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Abstract 



A model of quark mass matrices from six dimensions, which is nearly 
democatic in nature and which is previously constructed by two of us (PQH 
and MS), is studied in detail in this manuscript. We found that not only it fits 
all the six quark masses as well as the CKM matrix but also that there exists 
a region in the allowed parameter space of the model where the constraint on 
the parameter 9 of the Strong CP problem is satisfied. This region itself puts 
a constraint on the CKM parameters p and f]. As such, through our analysis, 
there appears to be a deep connection between Strong and Weak CP in this 
model. 
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I. INTRODUCTION 



The search for a plausible model of fermion masses is a continuing quest in particle 
physics. In particular, the quark sector is a fertile ground to test various models since it 
is there that one has the largest amount of information: quark masses, CKM angles, CP 
phase, and it is in this sector that most models turn their attention to. 

Two of us have recently constructed a model of fermion masses [1] in which the mass 
matrix is almost of the pure phase form and is constructed out of four plus two extra 
compact spatial dimensions. As shown in [1], one extra compact spatial dimension was 
needed to give a democratic mass matrix and another one was needed to make its matrix 
elements complex. In [1], an almost pure phase mass matrix was found to take the form 
M = gyv /\/2{l — Qjj)exp{idij) with i,j = 1,2,3, Qij « 1, and dij « 1, although Rcf. [1] 
contains a more general result. Our motivation for that work was based on an attempt to 
give a theoretical basis for the so-called pure phase mass matrix (PPMM) ansatz (similar 
to the previous form but with g^j — 0) of Ref. [2,3] which, at the time of its construction, 
was quite successful in fitting the various angles and masses. In so doing, we arrived at 
a mass matrix which contains the pure phase form as a particular limit. As we shall see 
below, the general result of Ref. [1] allows us to be able to fit the latest determination of 
the CKM elements [4]. Along the way, as stipulated in Ref. [1], we found a special region, 
in the allowed parameter space that fits the CKM matrix, where the parameter 9 of the 
famous strong CP problem can be found to satisfy the experimental bound 9 < 10~^. This 
result is somewhat surprising since it is not at all evident that solutions of our model that fit 
the mass spectrum and the CKM matrix could also give values of 9 below the experimental 
bound. This connection between weak and strong CP is certainly very intriguing and will 
be the subject of our focus at the end of the paper. 

The organization of the paper will be as follows. We first briefly review the construction 
of a democratic mass matrix (DMM) in five dimensions. We then summarize the salient 
points of the model of Ref. [1]: its construction in six dimensions and the resulting quark 
mass matrices. Next we enumerate and describe the parameters used in the numerical 
analysis of the mass matrices which is carried out in the section that follows. There we 
will show the allowed region in our parameter space where solutions are found to fit both 
the mass spectrum and the CKM matrix. Finally, we discuss a subspace in the allowed 
region where the bound 9 < 10^^ is obeyed. In particular, wc present some thoughts on the 
possible physics which might be responsible for this behaviour. 

II. DEMOCRATIC MASS MATRICES FROM FIVE DIMENSIONS 

Before discussing the results of Ref. [1], let us first review how a democratic mass ma- 
trix (DMM) [5] arises in the case with one extra compact dimension. A DMM is a spe- 
cial case of an almost PPMM with Qij = 0, 9ij = 0, namely one in which all matrix 
elements are unity, apart from a common factor, and hence the name "democratic". As 
described in [1], in order to obtain a democratic mass matrix and to avoid unwanted Flavour- 
Changing Neutral Current processes (FCNC), we imposed the following permutation sym- 
metries on the Action: ® Sf ® Sf, with Q S^Q, W S^^U" and S-^^D'. 
Q{x",y), [/'^(x",y), D'^{x",y) are the five-dimensional Dirac fields whose left-handed zero 



2 



modes are given respectively by q'(a:"), and d'j^{x"). (For convenience, left-handed 

fields are used throughout [1] and in this paper with the last two fields representing actually 
the two quark SU{2)l singlets.) The extra dimension is compactified on an S1/Z2 orbifold. 
The action which obeys this permutation symmetry is the sum of two terms, Sq and Syuki 
where 

So = j d'>xQ.,{iI?^ + f<P{y))Q, + Ut{iI?^ + fct>{y)-mu)Ut 

+Dt{iI^, + fcl>{y)-mn)Dt, (1) 



Svuk = / d'xku E QJC^H HU^+I d'xko E QJ^^H Y1 D'j + h.c. . (2) 

i j i j 

In Eq. (1) is the covariant derivative. (The gauge fields are supposed to be uniformly 

spread along the fifth dimension y inside the thick brane.) (f)(ii)5ij is the vacuum expectation 
value VEV for the background scalar field $(a;", y). The attractive proposal of [7] to localize 
chiral zero modes at different points along the extra dimension y was adopted in [1]. As 
a result, mu and mo are the five dimensional "masses" which determine the locations of 
ii^(x") and d'y^x"') along y. (As pointed out in Ref. [1], in order to have an invariant "mass 
term" under the Z2 symmetry, one has to require a "mass reversal" m —m. The behaviour 
of m under Z2 could come for example in a model where the "masses" are generated by the 
radiative corrections to the VEV 0(|/), (piy) — ^ (piy) + 50 with 54> being independant of 
y. Because $(a;",|/) — > — $(a;", L5 — 7/) under Z2 symmetry, at the same time (50 — > —50 
originating in this "mass reversal".) In Eq. (2) ku and ku are the Yukawa couplings in 
five dimensions which have been chosen real and flavor independant and H[x°',y) is the 
five dimensional SM doublet Higgs field whose zero mode h{x°') is assumed to be uniformly 
spread along y inside the thick brane. Here C5 = 707272/ is the charge conjugation operator 
in five dimensions. 

For the purpose of keeping track of the dimensionalities of various objects, let us remind 
ourselves that, in five dimensions, the Yukawa coupling ku has a (mass) dimension M~^/^. A 
scalar field, in five dimensions, has a dimension M^/^. The zero mode of the SM Higgs field 
can be written as K h{x^) where h{x°') is the usual 4-dimensional Higgs field with dimension 
M, and therefore the constant C has a dimension M^^"^. In dimensionally reducing the above 
action to four dimensions, the following dimensionless combination appears in the Yukawa 
term: 

9y,u = kuK , (3) 

The effective action for the Yukawa term of the Up sector in four dimensions can now 
be written as 

Syii = / d'xgY,u E qI{x-)h(x-)Cu%x-) J dyi\{y)il{y) + h.c. . (4) 

where a similar expression holds for the Down sector. As stressed in Ref. [1], since all the 
qis are located at the same place inside the brane, and similarly for all the ufs. The wave 
function overlap / dy^]^{y)^-lc{y) is independant of i, j. The Yukawa action now becomes 
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/I 1 1 

S^fJ^ = / d^xg^^^lq^{x")h{x") 1 1 1 I Cu^x'^) + h.c. 
' \1 1 1 

where gy^l is given by 

9y,1 = 9y,u j dy^q{y)^u-{y) ■ 
and similarly for the down sector 

QyJ = 9Y,d j dy^q{y)^d-{y) ■ 
From Eq. (5) one obtains the democratic mass matrix (DMM) 

/1 1 1\ 



1 1 1 

VI 1 ly 



(5) 



(6) 



(7) 



(8) 



which has eigenvalues 3gYv/V2,0,0, with v ~ 2A6GeV. The DMM of Eq. (8) does not 
reproduce the right mass spectrum and the right CKM matrix. What has been shown in 
Ref. [1] is that by adding another compact extra dimension one can obtain a viable scenario 
represented by an almost-PPMM. 



III. ALMOST PURE PHASE MASS MATRICES FROM SIX DIMENSIONS 

A. The Action 

The main idea of Rcf. [1] is that by introducing a sixth compact extra dimension, and 
by requiring that the background scalar field couples to the fermions through a Yukawa 
interaction which is non-local along that extra dimension, one can obtain an oscillatory 
behavior for the fermion wave function along the sixth dimension. Fermions are delocalized 
along the sixth dimension, in contrast with the fifth dimension, and the oscillatory behavior 
of the wave functions, together with the breaking of family symmetry, has the effect of 
producing phases in the mass matrix. Let us now summarize the main results obtained in 
Ref. [1]. We first rewrite more compactly the action given in Ref. [1] which is the sum of 5*0 
and Syuk where 

So^l (fx \Q,iI?^Q, + Q,{z){Ll>,{z) - ^)Qi{L, - z) - Q,{L^ - z){^-(t>i{z) - '^)Qi{z) 

+DtiI^,Dt + Dt{z){^-cl>,{z) - '^)Dt{L, -z)- Dt{L, - z){^-m - '^)Dt{z) 
+fQ^V,<P\y)Q, + UtT,{f<P\y) - mu)U^ + DlT,{f4>\y) - mD)Dl] , (9) 

and 
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(10) 



where Cg = ror2r2 is the charge conjugation in six dimensions. (The gamma matrices in 
six dimensions can be obtained in Ref. [1].) In Eq. (9) we expressed the dependance from 
z only for the non-local interaction terms. The important point here is that while these 
interactions will produce an oscillatory behavior for the fermion wave functions along the 
sixth dimensions, the local terms, which are built using Fy, are found to be responsible for 
localizing the fermions along the fifth dimension. The above actions are invariant under the 
family permutation symmetry (8) S^" (8) S^'^ . 

The vacuum expectation values (VEV) for the background scalar fields $ and are 
given by 



< ^{x'',y,z) >-- 




(11) 



and 



< ^'{x'',y,z) >-- 



( <P'{y) 

(p'iy) 
V 0'(y) 



(12) 



As in Ref. [1], the family symmetry is broken by the background scalar field (f)i{z)Sij and 
by introducing different non-local "mass terms" mQ^i, mjj^i and m/j j for each family. To 
break the family symmetry together with the left right symmetry along the sixth dimension 
will allow us to reproduce the right mass spectrum and the right CKM matrix. 

As shown in Ref. [1], the absence in Eq. (9) of local interactions of the form 'I'^^, 
which will localize the fermions wave function along the sixth dimension, is obtained by 
introducing the discrete symmetry Q 



■0(x", z) 

'^{x'^.Lq - z) 

$(x",z) 



Qi^{x'',z)^T,'4^{x'^,z), 

Q^ix'^, Le-z)^ -r^iPix'^, Le-z), 

Q^x'^,z) = $(x",^). 



(13) 



As pointed out in Ref. [1] the realization of the Q-symmetry of Eq. (13) implies that the 
introduced orbifold for the compactification is S\l{Zi x Zg) instead of Si/Z2. This also 
implies that the physical space is [0,L6/2] instead of the initial support [O.Le]. The non- 
local terms of Eq. (9) and the local terms containing Fy are invariant under the above 
Q-symmetry. 

B. The Mass Matrix 



Prom the Yukawa action of Eq. (10) one can now obtain the effective action in four 
dimensions 

Sy^i- j d'xgY,uj:<lIi^")hix")Cu%x") J dyil^{y)iiAy) j dzC,:,{z)iiA^) + h.c. . (14) 
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where we considered only the up sector, but equal considerations hold for the down sector. 
Using Eqs. (5) and (6) one can rewrite SyJ), as 



where 



( 



Oil ^12 Ol3 



O2I Q22 ^23 
\ 031 032 033 / 



1 fLe 

— — dzexp[—i{2fViln{cosh{iXiz))/iXi — 2fVjln{cosh{iXjz))/ijij 
Lq Jo 

-{mq,i-muj)z)] 
In Eq. (16) we have used for and ^g „c respectively 

1 



--exp[i{2fViln{cosh{niz))/ni - ruq^iz)] 



(15) 



(16) 



(17) 



ilu- = -^exp[i{2fViln{cosh{iiiz))/iii - mu,iz)] (18) 

V-^6 

which correspond to a VEV (pi{z) = Vjtgh^fiiz) with fii = ^\/2Vi. 

As pointed out in Ref. [1], Lq is now a generic symbol for the length of the physical space, 
which is Lq/2 for the orbifold Si/{Z2 x Z2). From Eq. (14) one obtains the mass matrix 

^ / ail ^12 ^13 
''^^ 021 023 I ■ (19) 



\ ^31 032 033 



Following Ref. [1], if one now uses the linear approximation for the kink, which is valid for 
l/ytij ~ 0{Lq), all domain wall thiknesses along z are of the size of the compact dimension, 
one can obtain for the elements a^j the form (1 — Qij)e^^^^ with g^j <^ 1 and 6ij -C 1. In the 
linear approximation for the kink one obtains the following expressions for 



1 f^e 



dij — J. 



' dz exp[-i{Anlz'^ - Aruij)] , (20) 



where 

A/.J = l(2/y,/.,-2/y,/.,), (21) 

Arriij = ruq^i - rriuj ■ (22) 
As shown in Ref. [1], one can explicitly carry out the integration and obtain 
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Now if ^J'Ajj^jLq = Xij <^ 1 and Arriij =<^ 1 one can expand Eq. (23) giving 

= (i - ^4 - ^4- + 1^4?/«.) e^-p (f - t)) + ' ^^^^ 

which has the desired almost-PPMM form. It has to be stressed here that the expression for 
the mass matrix which has been used to make our numerical simulation is the one from Eq. 
(16). This implies that when wc looked for a solution in the parameter space, we did not 
have to restrict ourselves to the particular range of values for the parameters where both 
the linear approximation for the kink and Eq. (24) were valid. 

As pointed out in Ref. [1], by looking at Eqs. (19) and (16) one can make the impor- 
tant following consideration. If one introduces the same "mass" term for left and right 
components for each family i, which means that rriq^i = rriu^i (and similarly mq^i = rrid^i for 
the Down sector), then the mass matrix of Eq. (19) is hermitian, i.e. aij = a*^. In this 
particular case the mass matrices for the Up and Down sectors differ only by the Yukawa 
couplings and one will not be able to reproduce a realistic mass spectrum. Moreover the 
diagonalization matrices are the same, i.e. Vu = Vu, and Vckm = ^J^d becomes just the 
unit matrix. So in order to obtain a realistic mass spectrum and CKM matrix one needs 
to introduce different "mass" terms for left and right components at least for one sector. 
What will be shown in the following is that we will be able to reproduce the right mass 
spectrum and the right CKM matrix for the case in which both up and down mass matrices 
are not hermitian. While hermitian mass matrices do not give a realistic scenario, they have 
the important property of having a real determinant. This implies that the argument of 
their determinant is zero. As pointed out in [1], this fact could form the seed for a solution 
to the strong CP problem. 



(23) 



IV. DESCRIPTION OF THE PARAMETER SPACE 

In this section we are going to describe the parameter space for the model of Ref. [1]. 
The particular case we consider has 10 parameters. What has to be said here is that we 
started our numerical simulation considering cases with a higher number of parameters and 
only after examining the results obtained, we were able to reduce the parameter space to 10. 
This, by the way, is the same number of parameters found in the quark sector: Six masses 
and four CKM parameters. The 10 parameters we are considering are the following: 

• gyl and gyl defined respectively by Eq. (6) and (7) 

• /ii, /i2 and whose inverses give the thickness of the domain walls 4>i{z), 4>2iz) and 
Mz) of Eq. (11) 

• Eq^i = rUq^i^ \/2/2f and 5^,2 = "^g,2\/A/2/2/ with ruq^i the "mass" term for the 1** 
family left component and mq^2 for the 2'^'^ and 3'''^ family left components 
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• £u,2 = ?7i„,2'\/ A/2/2/ and eu,3 — mu,3^JX/2/2f the "mass" term respectively for the 2"^^ 
and 3^*^ family right components, e^i — 



• Ae being the common split in my A/2/2/ of £d,i = £u,i — Ae with respect to eu,i 

As can be seen from the particular choice of parameters the background scalar potential 
does not break the left-right symmetry 0^ = = (jyf but on the other hand breaks the 
family symmetry 0j ^ (pj. The left-right symmetry is broken by choosing Eq^i ^ Su4,i and by 
choosing different Si for different indices i one breaks additionally the family symmetry. The 
choice £g,3 = £q,2 comes from the analysis of the same model in the case of 11 parameters, 
where all the Sq^i were different, and which gave as a result Sq^^ ~ £g,2- 



V. RESULTS FOR MASS MATRICES FROM SIX DIMENSIONS 

In this section we will present the results obtained for the parameter space and for the 
quantities of Table 1. 

We should point out that each solution that we found corresponds to a point in the 
parameter space with all the fitted quantities of Table 1 being in the experimental range 
[4,8]. In Figs. 1,2 and 3 we give the masses of the 6 quarks in GeV evaluated at the 
Mz scale for three different cases corresponding to three different ranges of arg{detM) = 
arg{detMu) + arg{detMd). It will be clear in the next section why it is interesting to look at 
the quantity arg{detM) when we present a scenario for a possible solution to the Strong CP 
problem. We choose to evaluate the running masses mq{fj,) at the scale /i — Mz, because 
the CKM matrix parameters V^^^^ are given at /j, — Mz- This is a common approach for 
quark mass matrix phenomenology. See Ref. [9] for a review of the running masses and the 
renormalization group equation that describes the evolution of the running quark masses 
mq{n) with the scale The edges of each box in Figs. 1-3 give the uncertanties for the 
masses, which depend not only on the errors of the input parameters for the renormalization 
group (RG) equation, but also on the error of the parameter which governs the flow itself, 
i.e. the strong coupling as{Mz). 

In the following we present two numerical examples corresponding respectively to 
arg{detM) ~ 0(1) and arg{detM) < 2 x 10~^°. For each numerical examples we will 
also give the corresponding parameter space. Below we rewrite the expression for the mass 
matrix to make clear the role of each parameter. 



9Yeff—^ — exp (-^2/y'^J ^ dz exp[-i{ln{cosh{iiiz)) - ln{cosh{iijz)) - {sq^i - eu,j)z) 



^/2. 
V 1 
72 ^ 



where £q^u,i = i^q,u,i\l ^/^/^f with /ij = ^\/2Vi. In our numerical simulation we will ignore 

the phase factor exp{—i 2f^J^) in Eq. 25, which is independent of the indices All results 
correspond to Le = 1. 



(25) 
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One word of caution is in order here. In comparing our results with the phenomenological 
extractions of the CKM matrix elements, we take into account the following points. 1) 
The magnitudes of Vub and Vcb are obtained from tree-level decays and are to a very good 
approximation independent of contributions from new physics. 2) If one were to use the 
Unitarity Triangle parameters p and fj for comparison, one has to assume that possible new 
physics contributions (through loop effects for example) conspire to bring the apex (p, fj) of 
the Unitarity Triangle to within the allowed band of the "unitarity clock" . In the following 
we will use both points 1 and 2 to make our comparisons with experiments. 

First, we give a numerical example corresponding to the case arg{detM) ~ 0(1). 
In Eq.s (26)- (30) we give the parameters space for the first case. 



= 7.378 , ii2 = 8.460 , 113 = 8.531 , 



(26) 



.262, £„2 = 5.090, £„3 = 5.090, 



£„i = 0.000, £„2 = 1.120, £„3= 1.198, 



(27) 
(28) 



£a = -2.044, £d2 = -0.924, £^3 = -0.846 , 



(29) 



gYuv/V2 = 152.31 , gYdv/V2 = 24.48 . 



(30) 



We have decided to present the parameter space in a more readable way but it is important 
to remember that the number of independent parameters is 10. In Eqs. (31) and (33) we 
give the numerical expressions for the up and down quark mass matrix for the first case and 
in Eq.'s (32) and (34) the corresponding mass eigenvalues. As one can see, the elements of 
the mass matrices, although not equal to each other in magnitudes, are of the same order, 
revealing their democratic origins. 



M, = 152.31 GeV 



/ 0.1111 + 0.1690i -0.1937 - 0.4167i -0.1852 - 0.4377i \ 

0.1039 + 0.1750i -0.1857 - 0.4221^ -0.1765 - 0.4430i 
V 0.1031 + 0.1758i -0.1846 - 0.4238i -0.1753 - 0.4446i / 



(31) 



m„ = 0.0023 GeV , = 0.624 GeV , rrit = 183.0 GeV , 



(32) 



Mrf = 24.48 Gey 



/ -0.0105 + 0.0003i -0.0424 - 0.0086i -0.0544 - 0.0121i \ 

-0.0090 + 0.0041i -0.0442 - O.OOGOi -0.0563 - 0.0099i 
V -0.0081 + 0.0043i -0.0451 - 0.0060i -0.0572 - O.OlOOi / 



(33) 



rud ^O.OO^S GeV , m, = 0.081 Gel/ , rub^ 3.09 GeV. 



(34) 



In Eqs. (35) and (36) we also give the absolute values of the mass matrices of Eqs. (31) 
and (33) which show that breaking the family symmetry does not destroy completely 
the democratic structure of the mass matrices. 
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Mu = 152.31 GeV 



/ 0.2023 0.4595 0.4753 \ 

0.2035 0.4611 0.4768 , 
V 0.2038 0.4622 0.4780 / 



(35) 



Mrf = 24.48 GeV 



0.0105 0.0432 0.0558 
0.0099 0.0446 0.0571 
0.0092 0.0455 0.0581 



(36) 



In Eqs. (37) and (38) we give the absolute values of the mass matrices of Eqs. (31) and 
(33) with rescaled values of the matrix elements. Eqs. (37) and (38) show in a more explicit 
way that the deviations from a democratic mass matrix are of 0(1). 



M„ = 72.80 061/ 



/ 0.4232 0.9614 0.9943 \ 

0.4257 0.9647 0.9976 , 
V 0.4264 0.9671 1.0000 / 



(37) 



= 1.42 GeV 



0.1807 0.7446 0.9601 
0.1702 0.7675 0.9836 
0.1584 0.7836 1.0000 



(38) 



In Eq. (39) we give the CKM matrix corresponding to the mass matrices of Eqs. 31 and 
(33), in Eq. (40) its absolute value, and in Eqs. (41) and (42) the values for the parameters 
p and fj, and for sin 2/3 and 7. 



0.9711 - 0.0884i 0.0930 - 0.2014i -0.0036 + O.OOlOi 
VcKM^ I -0.1185 -0.1872i 0.9734 - 0.0392i -0.0144 + 0.0398^ 
0.0095 - 0.0064i 0.0150 + 0.0380z 0.9986 + 0.0302i 



(39) 



VcKM 



/ 0.9751 0.2218 0.0038 
0.2215 0.9742 0.0423 
V 0.0423 0.0409 0.9991 



(40) 



p = 0.18, 77 = 0.35 



sin 2/3 = 0.72, 7 = 63.2°. 



with p and rj being defined as 



p = Re{VuaV:,V:,V^)/\V^V^\' , 



fj = Im{Vu,V:,V:,Vc,)/\VeaV:,\' . 



(41) 
(42) 

(43) 
(44) 



and sin 2/3 and 7 as 
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sin 2^ 



2r7(l-p) 



(45) 



^^n 90° . 1 / 2pf] ^ 
7 = 90° sin-^ ' \ 



(46) 



In Eq. (47) we give the values for the arg{detMu), arg{detMd) and for their sum. 

arg{detMu) = -1.5692 , arg{detMd) = 1.8643 , arg{detM) = 0.2951 . (47) 

The following numerical results corresponding to a value of arg{detM) < 10~^ are pre- 
sented in the same way as the above example. 



fxi = 7.365 , //2 = 8.456 , //g = 8.532 , 



(48) 



£51 = -8.189 , eg2 = 5.026 , £^3 = 5.026 , 



= 0.000, £„2 = 1.105, £„3 = 1-188, 



= -2.274, £d2 = -1.168, £d3 = -1.085, 



(49) 
(50) 
(51) 



gYuv/V2 = 144.63 , gYdv/V2 = 23.60 . 



(52) 



0.1154 + 0.1622i -0.1846 - 0.4383i -0.1740 - 0.4608i \ 
M„ = 144.63GeV^ | 0.1096 + 0.1662i -0.179 - 0.4367i -0.1683 - 0.459H , 
0.1090 + 0.1670i -0.1781 - 0.4382i -0.1671 - 0.4606i / 



(53) 



m„ = 0.0028 GeV , = 0.621 GeV , rrit = 177.9 GeV , 



(54) 



Md = 23.60 Gey 



/ -0.0609 + 0.0113i -0.0167 - 0.0046i -0.0294 - 0.0063i \ 

-0.0630 + 0.01480i -0.0144 - 0.0006i -0.0272 - 0.0025i 
\ -0.0622 + 0.0147i -0.0152 - 0.0005i -0.0280 - 0.0024i / 



, (55) 



rud = 0.0047 GeV , m, = 0. 105 GeV , rub = 2.9 GeV . 



(56) 



Mu = 144.63 Gey 



/ 0.1990 0.4756 0.4925 
0.1991 0.4720 0.4889 
\ 0.1994 0.4730 0.4899 



(57) 



M<i = 23.60 Gey 



/ 0.0619 0.0173 0.0301 \ 

0.0647 0.0144 0.0273 
V 0.0639 0.0152 0.0281 / 



(58) 
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0.4042 0.9656 1.000 \ 
M„ = 71.23 Gel/ | 0.4043 0.9583 0.9927 , (59) 
0.4049 0.9604 0.9948 / 



0.9568 0.2676 0.4655 \ 
Md = 1.53 Gey I 1.0000 0.2229 0.4223 . (60) 
0.9873 0.2350 0.4346 / 



-0.9681 + 0.1152i -0.1707 + 0.1425i 0.0038 + O.OOlli 
VcKM= I -0.1879- 0.1183 0.9741 - 0.01456i -0.0078 + 0.0384i | , (61) 
0.0068 - 0.0096i 0.0082 + 0.0367i 0.9991 + 0.0126i 



CRM — 



/ 0.9750 0.2223 0.0040 \ 

0.2220 0.9743 0.0392 
V 0.0118 0.0376 0.9992 y 



(62) 



p = 0.31 , 77 = 0.32 , (63) 

sin 2/3 = 0.77, 7 = 46.3°, (64) 

arg{detM^) = -1.55845421528 , arg{detMd) = 1.55845421485 , arg{detM) = 4.3 x 10"^° . 

(65) 

In Figs. 4, 5 and 6 we give the absolute values of the CKM matrix elements for the three 
different cases corresponding to three different ranges of arg{detM). The uncertainties for 
each element are given by the edges of the corresponding window. Fig. 8 shows the solutions 
for p and fj. The sharp edges for the solution patches are due to the constraints imposed 
on p and fj. Fig. 9 shows instead the solutions for sin 2/3 and 7. It has to be said here that 
the solutions appear in patches because of the way the minimization process works. One 
obviously can not exclude other solution patches. For example, could we try to make the 
minimization process follow different paths by toying with the input parameters, temperature 
and number of iterations, (see appendix A), and by changing the initial conditions. 
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TABLES 



Xi <Xi> - x7^^"|/2 



rriu 


2.33 X 10-3 


0.45 X 10-' 


rric 


0.685 


0.061 


mt 


181 


13 


rud 


4.69 X 10-3 


0.66 X 10-' 


rus 


0.0934 


0.0130 


rrib 


3.00 


0.11 


rriu/md 


0.497 


0.119 


ms/rud 


19.9 


3.9 


\Vud\ 


0.97485 


0.00075 


1 Vus 1 


0.2225 


0.0035 


\Vub\ 


0.00365 


0.00115 


\Vcd\ 


0.2225 


0.0035 


\Vcs\ 


0.9740 


0.0008 


\Vcb\ 


0.041 


0.003 


\Vtd\ 


0.009 


0.005 




0.0405 


0.0035 




0.99915 


0.00015 


p 


0.22 


0.10 


fj 


0.35 


0.05 


sin 2/3 


0.78 


0.08 


7 


59° 


13° 



TABLE I. Central values and uncertainties for the masses of the 6 quarks evaluated at Mz, for 
the two ratios mu/md and ms/md, for the absolute values of the CKM matrix elements, the CP 
parameteres p, fj, sin 2/3 and 7. 
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The parameter space with the 10 parameters described in the previous section and which 
corresponds to the three different cases for the the three different ranges of arg{detM), are 
given in Figs. 10, 11 and 12. It has to be said here that the solutions presented correspond 
to a two step procedure. First we find solutions by requiring /ii, 112 and /is to be larger than 
unity, because one can use in the physical space [0, Lg] the kink solution for (t)i{z) instead 
of the kink-antikink approximate solution Vtgh{fiz)tgh{fi{LQ — z)) [6] (the kink-antikink 
solution has the property of vanishing at both orbifold fixed point z = and z = Lq as 
required by the imposed boundary conditions to compactify on an S1/Z2 orbifold). Second, 
using initial conditions from the parameter space already found, we looked for solutions 
whic correspond to very small ranges of //i, ^2 and The reason we did this is because 
we noticed that the parameters which were more relevant to fit the quantities of Table 1 
were the "mass" terms of left and right components, i.e Eq^i, eu,i and ed,i respectively. As a 
result, in order to understand the dependence of the found solutions on the "mass" terms, we 
decided to restrict the range of /ii, 112 and The way one controls a range for a parameter 
consists simply in adding that parameter to the quantities one wants to fit, modifying the 
function / of Eq. (69) given in appendix A. For our numerical simulation we set Lg = 1. 
This implies that the introduced parameters /i and e have to be multipied by Lq for a general 
case. 

We have presented in this section an analysis of the quark mass matrices obtained in Ref. 
[1]. In this analysis, we have used 10 parameters and were successful in fitting all 6 quark 
masses and all the parameters of the CKM matrix. This is shown in Figs. (1-10) for three 
separate ranges of the quantity \arg{detM)\ which was defined above. In these figures, each 
dot in the scatter plot represents a solution which fits the masses and the CKM matrix. As 
one can see, the solutions which correspond to \arg{detM)\ ~ 10^^'^ and which will have a 
significance to the Strong CP problem, are within la and 1.5a of the so-called R-fit and 
Bayesian fit respectively, as presented in Ref. [4] . 



VI. POSSIBLE CONNECTION BETWEEN STRONG AND WEAK CP 

In this section, we will discuss a possible connection between the region of parameter 
space where arg{detM) < 10~^° and the Strong CP problem. 

It is well-known that the non-trivial vacuum of QCD generates a P and CP violating 
term in the Lagrangian of the form 

^0 = ^G.^G^'' (66) 

where 

= Oqcd + arg{detM) . (67) 

This term (66) gives a contribution to the electric dipole moment of the neutron [10], [11], 
with the current experimental limit [13] being 9 < 2 x 10~^°. The mystery of why 6 should 
be so small constitutes the so-called Strong CP problem. 

The most famous and elegant solution to the Strong CP problem is the Peccei-Quinn 
mechanism [14] where 6 becomes a dynamical field and relaxes to zero at the minimum of 
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its potential. This dynamical field manifests itself as a pscudoscalar particle- the so-called 
Axion- whose decay constant is now severely constrained by astrophysical and cosmological 
arguments [15]. Another solution involving a massless up quark [16] is largely disfavoured 
by studies of chiral perturbation theory. A third type of solution to the Strong CP problem 
which has no axion, is the Nelson-Barr type of mechanism [17] which assumes exact CP at 
tree level and whose mass matrices have real determinants, and, as a result, the strong CP 
problem does not exist. It can arise at loop levels and can be "under control". However a 
realistic model of this type is yet to be constructed. 

We have already mentioned that when we introduce the same "mass" term for left and 
right components (Eqs. (19) and (16)), i.e. m^^, = m^^, = 1^(1,1, the up and down quark 
mass matrices are hermitian and consequently their determinant is real [1]. arg{detMy) and 
arg{detMd) are separately zero. This situation suggests that the symmetry of the Lagrangian 
that makes arg{detM) = at tree level is the "left- right" symmetry of the components 
along the sixth dimension. If we also assume CP conservation at the Lagrangian level, this 
scenario would provide a solution to the strong CP problem. However this symmetry has 
to be broken because, as we have already mentioned, the case in which the quark mass 
matrices are hermitian does not reproduce the right mass spectrum and the right CKM 
matrix. One could imagine a scenario where the "left-right" symmetry is spontaneously 
broken. This will induce some loop corrections that will make the mass matrices deviate 
from hermiticity, reproducing the right mass spectrum and right CKM matrix. At the same 
time the breaking of this "left-right" symmetry could induce at loop levels a nonvanishing 
6. But one can envision a scenario were the deviation from hermiticity happens in such a 
way that arg{detMu) and arg{detMd), each being now of 0(1), can cancel each other so as 
to keep e <2x IQ-^^. 

As we have mentioned earlier, the solutions presented in this paper correspond to very 
small windows for the domain wall parameters //j's, because we wanted to put in evidence 
the effect of the "left-right" symmetry breaking, especially in the 6 parameter. In searching 
for a quantity which could "retain the memory" of this "left-right" symmetry, we decided to 
plot, see Fig. 13, the sum of the arguments \arg{det{M))\ versus the "CoM", the weighted 
average of the "mass" terms along the sixth dimension, which is defined as 

CoM = ^ ^I'i + ^^'i + .ggx 

12 ■ ^ ' 

As one can notice from Fig. 13, the sum of the arguments tends to go to zero for a particular 
CoM ~ 0.125, and this behavior is confirmed in Fig. 14. Notice that, for practical reasons, 
we have reduced the number of points (i.e. the number of solutions) in Figs. 13 and 14 
in order to reduce the sizes of the files containing these two figures. The actual number of 
solutions is much larger than what is shown in these figures. 

Now one can think that the value of 0.125 for the "CoM" when the "left-right" symmetry 
is broken, and which corresponds to ^ < 2 x 10"^'^, was also the value for the "CoM" before 
the breaking, when 6 was equal to zero. In other words one can imagine a scenario where 
the "mass" terms of left and right components are split such a way as to retain the same 
value of the "CoM". To invent a mechanism which could break the "left-right" symmetry 
and which could reproduce the scenario mentioned above is beyond the scope of this paper, 
but it is one of the main topics we would like to investigate in the future. 
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If one now looks at Fig. 8 which show the solutions in the p, fj plane for \arg{detM) \ in the 
three different ranges, one can see that the solutions corresponding to 10^^^ < \arg{detM)\ < 
2 X 10~^° tend to favor a particular region of the p, fj plane. In particular, the region of allowed 
solutions shrinks down in the p direction constraining the parameter p in the small window 
~ [0.3, 0.32]. These solutions are within 1 a and 1.5 a of the so-called R-fit and Bayesian fit 
respectively as can be seen from Fig. 8. This could be an artifact of the way the minimization 
process works, but it could also be an indication that there is a deep connection between 
Weak CP violation and Strong CP violation in our model. 

As we have already mentioned earlier, one cannot exclude completely the existence of 
other solutions, because of the way the minimization process works, but we believe that 
the results obtained here can give some very good indications of how the allowed parameter 
space might look like. 

VII. EPILOGUE 

We have presented in this paper a complete phenomenological analysis of the model of 
quark mass matrices as derived from six dimensions by [1]. With just 10 parameters, we 
have found a large number of sohitions which can fit the 6 quark masses as well as the 
CKM matrix as can be seen in Figs. (1-9). Furthermore, we subdivide these solutions 
into three sets: 1) Those that have 10^^ < \arg{detM)\ < 5 x lO^"*^, 2) those that have 
10"'^ < \arg{detM)\ < 10-^ and finally 3) those that have 10"^^ < \arg{detM)\ < 2 x 10"^°. 
The first two sets are given for the sole purpose of comparison with the last set which is most 
relevant to the Strong CP problem as discussed in the last section. As one can observe from 
Figs. (8) and (9), there is a deep connection between the weak CP parameters and the Strong 
CP phase. In order to satisfy the constraint on the 9 parameter, i.e. \arg{detA{)\ < 2 x 10^"*^° 
in our framework, the solutions obtained for p and f] are found to be within 1 a and 1.5 cr of 
the so-called R-fit and Bayesian fit respectively. As measurements of weak CP parameters 
become more accurate in the future, they will either rule out or confirm our predictions. 

ACKNOWLEDGMENTS 

We would like to thank Andrzej Buras and Gino Isidori for valuable comments. One 
of us (PQH) would like to thank the theory group of Laboratorio Nazionale di Frascati for 
hospitality during the course of this work. (AS) would like to thank Mike Timmins for 
introducing the simulated anneahng method and Ngoc-Khanh Tran for useful! discussions. 
Our work is supported in parts by the US Department of Energy under grant No. DE-A505- 
89ER40518. 

VIII. APPENDIX A 

As mentioned above the model we considered has 10 free parameters. Relying on this 
freedom we were able to fit the 6 quark masses evaluated at an energy scale equal to the 

mass of Z gauge boson M^, with constraints for the two ratios rriu/md and nis/md, the 
absolute values of the CKM matrix elements and the CP parameters p, fj, or equivalently 
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the three angles and one phase of the CKM matrix standard parametrization, for a total 
number of 10 quantities. The approach we used to derive the parameter space consists in 
minimizing a particular function, built in such a way that its global minima correspond to 
the region defined by the experimental constraints. This function is defined in the following 
way: 



where 9{x) is the step function, N is the number of quantities that we want to fit, xf^ is 
the predicted value for the ith quantity, x™" and .x™""^ fix the range for the ith quantity, 
and < Xi > is its average value. It is immediate to verify from Eq. (69) that when all 
the predicted quantities xf^^s are contained in the proper ranges, the function / takes its 
minimum value equal to zero. 

The set of parameters which correspond to a zero value for the function / is called a 
solution. In our particular case the function we are considering is a mapping from 7^^° to 7^. 
The parameter space is really big and to find the solutions, which correspond to the global 
minima of /, can be very challenging. The global minima of / can be infact surrounded 
by a lot of local minima which most of the time can make the minimization process fail. 
Therefore, just trying different initial conditions with the hope to find the good one that 
will allow the minimization process to follow the right path towards a global minimum can 
result quite inefficient. 

It is a common belief that there does not exist a general recipe to follow for minimization 
problems. A minimization procedure that can be very efficient for a particular problem, can 
be very inefficient in an another case or even fails. In the particular case where the function 
that we want to minimize depends on many parameters, there is a minimization procedure, 
called simulated annealing [18] [19], which seems to work more efficiently than others. This 
procedure is mostly used when the global minima are surrounded by a lot of local minima. 
Infact this minimization process can find a global minimum also after being trapped in a 
local minimum. The way instead most of the minimization processes work is to go, from the 
starting point, immediately downhill as far as they can go, but this often leads to a local 
minimum. 

At the heart of the method of simulated annealing is an analogy with thermodynamics, 
specifically with the way liquids crystallize, or metal cool and anneal. When a liquid is 
cooled down sufficiently slowly, the atoms are often able to line themeselves up and form 
a pure crystal, which corresponds to the state of minimum energy for the system. But if 
the liquid is cooled quickly, it does not reach this state but rather ends up in an amorphous 
state of higher energy. So nature, as long as the process of cooling is sufficiently slow, is 
able to find the minimum energy state. The way nature works is based on the fact that for 
a system in thermal equilibrium at temperature T, the probability for the system to be in 
a state of energy E is given by the Boltzmann probability distribution: 




< Xi> 



1=1 




< Xi> 



(69) 



i=l 




(70) 



17 



This implies that even at low temperature there is a chance, albeit a tiny one, for the 
system to be in a high energy state. Therefore, there exists the possibility for the system 
to get out of a local energy minimum and move towards the global one. This principle 
has been incorporated in what is called a Metropohs algorithm [20]. Given a simulated 
thermodynamic system, it is assigned a probability P — exp{—{E2 — Ei) / KT) to the change 
from a configuration with energy Ei to one with energy E2. If E2 < Ei, P is greater than 
unity and in this case to the change is assigned a probability P = 1, which is equivalent to 
say that the system always makes such a change. In the case E2 > Ei, one can compare the 
probability P — exp{—{E2 — Ei) / KT) with a random — number and make the change to the 
new configuration only if P > random — number. The system always takes a downhill step 
while sometimes takes an uphill step. The Metropolis algorithm can be used for systems 
other than the thermodynamic ones if we give: 

• A description of possible system configurations. 

• A generator of random changes in the configuration. 

• A function E (analog of energy) whose minimization is the goal of the procedure. 

• A control parameter T (analog of temperature) and an annealing schedule which tells 
how it is lowered, e.g. after how many random changes in configuration and with 
which step. 

Going back to our particular case, the function / we want to minimize is the analog of 
the energy, and each possible set of parameters correspond to a particular configuration of 
the system. For the algorithm to work a control parameter T, with an annealing schedule 
by which it is gradually reduced, has to be introduced, as well as a generator of random 
changes in the configuration, that is in the parameter space. The way these random changes 
are taken is the following: a positive quantity , given by — T • ln{random — number), is 
added to the stored function value. The same quantity is subtracted from the function value 
corresponding to every new set parameters that are tried as a replacement point in the 
parameter space. (The new points are obtained using the downhill simplex method. The 
simplex is a set of + 1 points with N the number of parameters, and the changes happen 
through reflections, expansions and contractions of the simplex) . As mentioned before, this 
method allowed the system to jump from a local minimum and to look for a global one. The 
algorithm that we have used has been taken from numerical recipes [21]. Other than the 
initial set of values for the parameters, also the value for the temperature and the number 
of iterations which determine the annealing schedule to reduce the temperature, has to be 
given as input. The way these last two values were chosen, as well as the initials conditions 
is mostly the result of different attempts. The output of our code is a patch of solutions, 
which have been subsequently tested using an independent code written in Mathematica. 
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FIGURES 



FIG. 1. Solutions for the 6 quark masses corresponding to 10""^ < \arg{detM)\ < 5 x 10~^ . 
The masses in GeV are evaluated at the Mz scale. The range for each mass as given in Table 1 is 
defined by the edges of the corresponding window. 

FIG. 2. Solutions for the 6 quark masses corresponding to 10^^ < \arg{detM)\ < 10^^ . The 
masses in GeV are evaluated at the Mz scale. The range for each mass as given in Table 1 is 
defined by the edges of the corresponding window. 

FIG. 3. Solutions for the 6 quark masses corresponding to 10~^^ < \arg{detM)\ < 2 x 10~^^ . 
The masses in GeV are evaluated at the Mz scale. The range for each mass as given in Table 1 is 
defined by the edges of the corresponding window. 

FIG. 4. Solutions for the absolute values of the CKM elements corresponding to 
IQ-^ < \arg{detM)\ < 5 x 10"^ The range for each element as given in Tabic 1 is defined by 
the edges of the corresponding window. The ranges delimeted by the dashed lines correspond to 
the new evaluations for \Vus\ (region on the left of the dashed line), \Vub\ and Vcb as in Ref. 4. 

FIG. 5. Solutions for the absolute values of the CKM elements corresponding to 
10~^ < \arg{detM)\ < 10"^. The range for each element as given in Table 1 is defined by the 
edges of the corresponding window. The ranges delimeted by the dashed lines correspond to the 
new evaluations for |T4s| (region on the left of the dashed line), \Vub\ and Vcb as in Ref. 4. 

FIG. 6. Solutions for the absolute values of the CKM elements corresponding to 
10^12 < \arg{detM)\ < 2 x 10"^°. The range for each element as given in Table 1 is defined 
by the edges of the corresponding window. The ranges delimeted by the dashed lines correspond 
to the new evaluations for \Vus\ (region on the left of the dashed line), \Vub\ and Vcb as in Ref. 4. 

FIG. 7. Solutions for the ratio |^6|/|V^| for the three cases corresponding to a) 
10"^ < \arg{detM)\ < 5 x 10'^, b) 10"^ < \arg{detM)\ < 10'^, and c) 
10-^2 < \arg{detM)\ < 2 x 10"^°. 

FIG. 8. Solutions for p and f] for the three cases, from top to bottom, correspond- 
ing to a) 10"^ < \arg{detM)\ < 5 x 10'^, b) 10"^ < \arg{detM)\ < 10~\ and c) 
10-^2 < \arg{detM)\ < 2 x lO"^'^. The delimeted area is the allowed region in the p, rj plane 
as obtained from Ref. 4. 

FIG. 9. Solutions for sin 2/3 and 7 for the three cases, from top to bottom, correspond- 
ing to a) 10"^ < \arg{detM)\ < 5 x lO'^, b) lO'^ < \arg{detM)\ < 10'^, and c) 
10-^2 < \arg(detM)\ < 2 x 10-^°. 
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FIG. 10. Summary of the 10 parameters space corresponding to 10 ^ < \arg{detM)\ < 5 x 10 ^. 
FIG. 11. Summary of the 10 parameters space corresponding to 10~^ < \arg{detM)\ < 10"^. 

FIG. 12. Summary of the 10 parameters space corresponding to 10~^^ < \arg{detM)\ < 
2 X 10-1°. 

FIG. 13. Solutions for \arg{detM)\ in the region [10-3, 5 X 10-1] vs the CoM as defined by Eq. 

68. 

FIG. 14. Solutions for \arg{detM)\ in the region [10-^2,2 x 10"!°] vs the CoM as defined by 
Eq. 68. 
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